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Abstract 

We have searched for the K + — > 71+77 decay in the kinematic region with ir + 
momentum close to the end point. No events were observed, and the 90% confidence- 
level upper limit on the partial branching ratio was obtained, B(K + — > 7r + 77, 
P > 213 MeV/c) < 8.3 x 10~ 9 under the assumption of chiral perturbation theory 
including next-to-leading order "unitarity" corrections. The same data were used 
to determine an upper limit on the K + — > 7r + 7 branching ratio of 2.3 x 10 -9 at the 
90% confidence level. 
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We report the results of a search for the rare decay K + — > 7i +/ ~f~f in the 
7r + momentum region P > 213 MeV/c from the E949 experiment [1] at the 
Alternating Gradient Synchrotron (AGS) of Brookhaven National Laboratory. 
The first observation of the decay in the ir + momentum region 100-180 MeV/c 
was reported [2] by the E787 experiment at the AGS with a partial branching 
ratio of B(K+ -> tt+77, 100 MeV/c < P <180 MeV/c) = (6.0 ± l.b(stat) ± 
0.7 (syst)) x 10~ 7 . In the region P >215 MeV/c no K + — > 7r + 77 decays were 
observed and, assuming a pure phase-space kinematic distribution, a 90% 
confidence-level (C.L.) upper limit of 5.0 x 10~ 7 was set on the total branching 
ratio [2]. This established that the K + — > 7r + 77 background to the rare decay 
K + — > 7T + z/z/ [3] was negligible. 

In an effective-field approach to low-energy hadronic interactions called chiral 
perturbation theory (ChPT) [4], there is no 0(p 2 ) ("tree-level") contribution 
to K + — > 7r + 77 or its neutral counterpart K\ — > 7r°77; the leading contribu- 
tions start at 0(p 4 ) [5]. For K + — > 7T + 77, both the branching ratio and the 7r + 
spectrum shape at 0(p 4 ) are sensitive to the undetermined coupling-constant 
c. There is no complete calculation at the next-to-leading order, 0(p 6 ). The 
dominant effects [6] are one-loop "unitarity" corrections, deduced from an 
empirical fit to the decay amplitude of K + — > 7r + 7r + 7r~ and containing the 
same constant c, and vector-meson exchange. In K + — > 7r + 77 vector-meson 
exchange is expected to be negligible compared to unitarity corrections. The 
corrections result in a slightly different prediction for the 7r + spectrum (Fig. 1). 
Separate fits to the measured n + spectrum from [2] with and without the uni- 
tarity corrections yielded c = 1.8 ± 0.6 and c = 1.6 ± 0.6, respectively, but 
slightly preferred their inclusion. For K° L — > 7r°77, the amplitude at 0(p A ) is 
determined without any undetermined coupling-constant, but the measured 
branching ratio, (1.41 ±0.12) x 10 -6 [7], is twice as large as predicted at 0(p A ); 
the vector-meson contribution in the next-to-leading order calculation (some- 
times parametrized by an effective coupling constant a v ) is considered to be 
important to this decay [8]. 

One of the consequences of the unitarity corrections to K + — > 7r + 77 is a 
non-zero amplitude in the kinematic region close to the end point of P =227 
MeV/c, where the two-photon invariant mass m 77 =0 MeV/c 2 , as shown in 
Fig. 1. The partial branching ratio B(K + — > 7r + 77, P > 213 MeV/c), cor- 
responding to m 77 <108 MeV/c 2 , is predicted to be 6.10ig;}2 x 10 -9 for c = 
1.8 ± 0.6 including unitarity corrections and 0.49^o; 2 | x 10~ 9 for c = 1.6 ± 0.6 
without the corrections. Observation of the decay at a partial branching ra- 
tio larger than predicted would indicate the contribution of vector-meson 
exchange or other new dynamics to K + — > 71+77. The kinematic region 
close to the end point in K® — > 7r°77 is known to be crucial to under- 
stand the CP-conserving component to the K\ — > 7r°e + e~ decay through 
the K\ — > 7r°77 — > 7r°e + e~ amplitude, but experimental results on a v in 
K\ — > 7T°77 [9,10] are inconsistent (see Ref. [11,12] for theoretical discus- 



3 



sions) . 



E949 was designed primarily to measure the decay K + — > -k + vu [13]. The AGS 
delivered kaons of 710 MeV/c to the experiment at a rate of 12.8 x 10 6 per 
2.2-s spill. Kaons, detected and identified by Cerenkov, tracking, and energy- 
loss counters, were slowed by BeO and active degraders, and came to rest and 
decayed in a scintillating-fiber target. The exposure of kaons in each spill was 
recorded by a scaler which counted the number of kaons entering the target 
during the time the data acquisition system was ready to accept data. Fig. 2 
shows a diagram of the apparatus. Measurements of charged decay products 
were made using the target, a central drift chamber, and a cylindrical range 
stack (RS) composed of 18 layers of plastic scintillator with two embedded 
layers of tracking chambers. The RS counters in the first layer (T-counters 
in Fig. 2) were 0.635-cm thick and 52-cm long; the subsequent RS counters 
were 1.905-cm cm thick and 182-cm long. The pion from the K + — > 7T + 77 
decay was identified by observation of the 7r + — > fi + — > e + decay sequence 
in the RS using 500-MHz waveform digitizers based on flash analogue-to- 
digital converters [14]. Trigger counters surrounding the target (I-counters in 
Fig. 2 which were 0.64-cm thick) and the T-counters surrounding the drift 
chamber defined the fiducial region. Counters (0.95-cm thick) surrounding the 
RS helped to suppressed the muons from K + — > and K + — > decays 
by identifying long-range muons that had completely traversed the RS. A 
hermetic calorimeter system surrounded the central region. The photons from 
K + — > 7r + 77 were detected in a lead/scintillator sandwich barrel detector 
(BV) surrounding the RS, while two endcap calorimeters and other subsystems 
(collar, "UPV", "AD", and "DPV" in Fig. 2) were used for detecting extra 
particles. A solenoid surrounding the BV provided a 1 T magnetic field along 
the beam line. 

The AGS proton-beam intensity, and the E949 beam line and apparatus were 
improved over those used in E787 [15] for the K + — > 7r + 77 study, which 
was performed in 1991 s . The new beam line [16] incorporated two stages of 
electrostatic particle separation which improved the acceptance for kaons as 
well as the K + to rr + ratio. The target, central drift chamber, and RS tracking 
chambers were replaced by a new target consisting of 0.5-cm square fibers, 
a new low-mass drift chamber [17], and straw-tube chambers, respectively. 
One third of the RS scintillation counters were replaced to increase the light 
output. A new photon detector, the barrel veto liner (BVL), was installed to 
add 2.3 radiation lengths of lead/scintillator sandwich material to the BV. The 
endcaps were replaced by new fully active detectors consisting of undoped-CsI 
crystals [18] with significantly increased light output, and both the target and 
endcaps were read out using 500-MHz waveform digitizers based on charged 
coupled devices [19] to improve timing and double-pulse resolution. Additional 



Many of the improvements were made in the E787 apparatus after 1991. 
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ancillary photon-veto systems [20] and a flasher system of Light Emitting 
Diode to aid in the RS energy calibration were also introduced. 

The 2002 data set used for this analysis derived from a total exposure of kaons 
entering the target N K = 1.19 x 10 12 . The trigger required a kaon at rest in 
the target to decay at least 1.5 ns later into a positively charged particle which 
subsequently came to rest in the RS, accompanied by coincident electromag- 
netic showers in both the BVL and BV, and no extra energy in the endcap or 
RS counters. The decay particle was required to penetrate the RS to at least 
the 16th layer to suppress backgrounds from the monochromatic K + — > 7r + 7r° 
decay (K n2 ) with P=205 MeV/c, and no further than the 17th layer to sup- 
press decays into muons. In the RS counter where the particle came to rest, 
called the "stopping counter" , a n + — > fi + decay was identified online based 
on the pulse shape information from the waveform digitizers on the RS. An 
improved trigger system [21] including a programmable trigger board reduced 
the online dead time. A total of 1.1 x 10 7 events met the trigger requirements. 

The momentum, equivalent range in plastic scintillator (R) , and kinetic energy 
(E) of the charged track were reconstructed with information from the target, 
drift chamber and RS. Tracks were accepted for the region defined by 213 
MeV/c < P < 234 MeV/c, 33.5 cm < R < 41.3 cm, and 116 MeV < E < 
135 MeV, where the lower limits corresponded to 3.3, 2.3, and 2.6 standard 
deviations, respectively, above the K W 2 peak (P = 205 MeV/c, R = 30.4 
cm, and E = 108 MeV). The larger search region (P > 213 MeV/c) compared 
with E787 (P > 215 MeV/c) resulted from improved kinematic reconstruction, 
which also removed the requirement of a constrained fit for consistency with 
K + — > 7T+77 kinematics that was used in E787 [2]. 

The timing and energy (E^) of the photons were determined by grouping 
adjacent hit modules in the BVL and BV to classify isolated photon showers 
("clusters"). The hit position in each module along the beam axis (z) was 
calculated from the end-to-end time and energy differences; the azimuthal 
angle (0) of the hit position in the end view of the detector was determined 
by the segmentation of the modules. The location of the photon shower in z 
and was obtained by an energy-weighted average of the hit positions and 
was used, in conjunction with the kaon-decay vertex position in the target, to 
determine the polar and azimuthal angles of the photon as well as the opening 
angles between the photon and the charged track in the side view (9 n +^) and 
in the end view In approximately half of the K + — > 7r + 77 decays the 

two photons were unresolved and appear as a single cluster in the BVL and 
BV, since the opening angle between two photons from K + — > 7r + 77 gets 
smaller for the events with tt + momentum close to the kinematic end point. 
The events with either one or two clusters were accepted; at least one cluster 
with 50 MeV < £ 7 < 320 MeV, w+7 > 155°, and ff+7 > 155° was required. 
For events with two clusters, the lower-energy cluster had to have E 7 > 10 
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MeV. 



Background sources from kaon decays at rest were classified into three types: 

"mismeasured"''' : K n2 decays with mismeasurements of the n + and the two 
photons, 

"of er/ap" : K n2 decays with the lower-energy photon overlapping the n + 

track in the RS, causing the kinetic energy of reconstructed 7r + track to 

be incorrectly measured, and 
"mwon" : K + decays with a yU + misidentified as a 7r + and with photons in 

the final state (e.g. K + — > K + — >• it n + v and fT^ with 7r + decay in 

flight). 

Another background source was due to kaon decay in flight: 

"DIF" : K n 2 decay-in-flight before the kaon came to rest in the target. 

Beam-related backgrounds (e.g. multiple beam particles scattering into the 
detector) were found to negligible. 

These backgrounds were studied from the data by imposing offline selection 
criteria ("cuts"). The requirements on the 7r + momentum, range and kinetic 
energy provided large suppression of all of the backgrounds from kaon decays 
at rest. Several cuts, referred to as "7 selection cuts" , were imposed to suppress 
the mismeasured background. These included cuts on the following quantities: 

• the cut on the invariant mass of two photons to reject events with two 
clusters and m 77 >100 MeV/c 2 , 

• the cut on additional coincident energy to reject events with activity not 
associated with the ir + and the candidate signal photons 9 , and 

• the cut on the photon clusters to reject events with two photons from a 7r° 
which hit the same or adjacent modules in BVL and BV and form a single 
cluster 10 [22]. 

In addition, only those events whose total photon-shower energy was deposited 
in the BVL and BV calorimeters were accepted. Cuts on dE/dx in the RS to 
reject events with a RS counter in which the measured energy was larger 
than expected from the reconstructed range in that counter were imposed to 
suppress the overlap background [22]. The cuts on the relation between the 

9 The extra activity was identified in the various subsystems, including the RS, 
as hits in the counters in coincidence with the tt + track within a few ns and with 
energy above a threshold of typically 1 MeV. 

10 Due to the kinematics of iVr2 and subsequent tt° — > 77 decays, the two photons 
must hit the modules at different z positions along the beam axis. An event was 
rejected if the maximum discrepancy among the z-position measurements in the 
modules of the cluster was larger than 113 cm. 
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range measured in the RS and the momentum measured in the drift chamber 
as well as the cuts on the ir + — > /i + — > e + decay sequence, recorded in the RS 
stopping counter, were also imposed to suppress the muon background. The 
DIF background was suppressed by requiring a delay of at least 2 ns between 
the time of the kaon coming to rest in the target and its subsequent decay, 
and by imposing the cuts on the timing between the 7r + in the RS and the 
K + in the Cerenkov counter. 

To study and measure these backgrounds, two independent sets of cuts were 
established for each background source. At least one of these cuts was inverted 
to enhance each background collected by the K + — > 71+77 trigger as well as to 
prevent candidate events from being examined before the background studies 
were completed. To avoid contamination from other background sources, all 
the cuts except for those being studied were imposed on the data. Possibilities 
of a correlation between the two sets of cuts or of a biased estimate of the 
effectiveness of the cuts were studied, and were found to be negligible. The 
signal acceptance and background levels were studied as a function of cut 
severity. A comparison of the observed background levels near but outside 
the signal region was made to the predicted background in these regions. All 
of these studies [23] were performed in the same manner as the K + — > 7i + uu 
analyses of E787 [24] and E949 [13]. Table 1 summarizes the background levels 
measured with the final analysis cuts and the two sets of cuts for studying each 
background source. In total, 0.197 ± 0.070 background events were expected 
in the signal region. 

The acceptance (A) and the single event sensitivity (SES) for K + — > 7r + 77 
in the kinematic region P > 213 MeV/c were derived from the acceptance 
factors in Table 2 and the total kaon exposure Nk = 1.19 x 10 12 times the 
fT + -stopping efficiency, which is the fraction of kaons entering the target that 
came to rest before decaying. The efficiency was measured to be 0.754 ±0.024 
with the K W 2 events collected with the K + — > 7r + 77 trigger and selected us- 
ing all the analysis cuts except those designed to remove K n 2 events. Thus 
the sensitivity for K + — > 7r + 77 was normalized to K v2 and many systematic 
uncertainties in the measurement of the acceptance factors, in particular in 
the trigger acceptance, canceled. We obtained A = (2.99 ± 0.07) x 10~ 4 and 
SES = (3.72 ± 0.14) x 10~ 9 for c — 1.8 including unitarity corrections and 
A = (1.10 ± 0.04) x 10~ 4 and SES = (1.01 ± 0.05) x 10~ 8 for c = 1.6 with- 
out the corrections. The former sensitivity was below the predicted branching 
ratio of 6.10 x 10~ 9 , giving an expectation of 1.6 events. In order to verify 
that the sensitivity estimations were correct, a sample of K + — > ji + u decays 
accumulated with a calibration trigger was analyzed. The branching ratio of 
0.628 ± 0.020 measured with the same fC + -stopping efficiency was consistent 
with the Particle Data Group value [7]. The systematic uncertainty in the 
sensitivity was estimated to be less than 10%. 
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After imposing all analysis cuts, no events were observed in the signal region 
(Fig. 3). The group of 74 events around R = 32 cm and E = 110 MeV are 
due to the K n2 background. Taking 2.24 events instead of zero according to 
the unified approach [25] with the background contribution of 0.197 events, 
we set a 90% C.L. upper limit on the partial branching ratio B(K + — > 7r + 77, 
P > 213 MeV/c) as 8.3 x 10~ 9 for c = 1.8 including unitarity corrections and 
2.3 x 10~ 8 for c — 1.6 without the corrections. The systematic uncertainty 
was not taken into consideration in deriving the limits. For the purpose of 
comparison with the previous E787 results, a 90% C.L. upper limit for the 
total K + — > 7r + 77 branching ratio assuming the phase-space distribution was 
calculated; the present limit 6.0 x 10" 8 is 8.3 times lower than the same limit 
in E787 (5.0 x 10~ 7 ). 

The data described above were also used to set an upper limit on the branching 
ratio for K + — > n + j decay, which is forbidden by angular-momentum conser- 
vation and gauge invariance, but is allowed in noncommutative theories [26]. 
The signature of K + — > 7r + 7 was a two-body decay of a kaon at rest with a 227- 
MeV/c 7r + track in the RS and a 227-MeV photon emitted directly opposite 
to it and observed in the BVL and BV calorimeters. The trigger, event recon- 
struction, and offline selection criteria in the study of K + — > 7r + 77 had been 
designed so that the same data were available to the search for K + — > tt + j. 
Since the background levels were already small, the 7r + accepted region was 
not reduced for the K + — > 7r + 7 analysis. The previous limit from the E787 
study was 3.6 x 10" 7 [22] (90% C.L.) from data collected in 1996-1997 with a 
highly prescaled trigger with relaxed trigger-conditions resulting in the total 
kaon exposure of 6.7 x 10 8 . The new 90% C.L. upper limit from E949, using 
the acceptance for K+ -> vr+7 of (1.08 ± 0.02) x 10" 3 , is 2.3 x 10" 9 . 

The results from this study cannot confirm nor rule out the unitarity correc- 
tions of ChPT, but the upper limits obtained are the the most restrictive yet 
achieved on K + — > 7r + 77 and as well as on K + — > 7r + 7. The E949 experiment 
has been shown to be suitable for the study of K + — > 7r + 77 in the n + momen- 
tum region close to the end point. The experimental uncertainty is limited by 
statistics; additional data is required to more stringently test the predictions 
of ChPT for this decay. The possibility to observe the K + — > 7r + 77 decay in 
the kinematic region, if the ChPT including unitarity corrections is correct, 
gives further impetus for additional data collection. 
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tt + momentum (MeV/c) 



Fig. 1. Predictions for the ir + momentum for the ChPT parameter c = 1.8 including 
unitarity corrections (solid line) and for c = 1.6 without the corrections (dashed line) 
as described in the text. The ratios of the partial branching ratio in the kinematic 
region P > 213 MeV/c, indicated by the arrow, to the total branching ratio are 
5.77 x 10~ 3 and 0.515 x 10~ 3 , for c = 1.8 and c = 1.6 respectively. 




Fig. 2. Schematic side (a) and end (b) views of the upper half of the E949 detector. 
C: Cerenkov counter; B4: energy-loss counters; I and T: inner and outer trigger 
scintillation counters; RSSC: RS straw-tube tracking chambers. New or upgraded 
subsystems for E949 (shaded) included the barrel veto liner (BVL), collar, upstream 
photon- veto (UPV), active degrader (AD), and downstream photon- veto (DPV). 
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Fig. 3. Range (equivalent cm of plastic scintillator) vs. kinetic energy (MeV) plot of 
the events with all analysis cuts imposed. The box indicates the signal region for the 
— > 7r + 77 decay. The dark points represent the data. The simulated distribution 
of expected events from K + — > tt + jj for c = 1.8 including unitarity corrections is 
indicated by the light dots. 
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Source 


background level 


two sets of cuts 


mismeasured 


0.017 ± 0.00b 


7r T accepted region (F,H,E) 


7 selection 


overlap 


O.Ooo ± O.Ooo 


7r + accepted region (P,R) 


7T+ dE/dx 


muon 


0.090 ± 0.020 


7r + accepted region (P,R,E) 
ir + range-momentum relation 




DIF 


0.025 ± 0.014 


delay in the target 


RS - C timing 



Table 1 

Expected background levels in the signal region and the two sets of cuts for studying 
each background source. The total background was 0.197 ± 0.070 events. 
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Acceptance factors UC w/o UC samples 



Trigger 


0.0623 


n n /i c\'~7 

0.0407 


ML, 


Kb 


7r T reconstruction and nducial cuts 


0.980 


n no r 

0.935 


ML, 




7r T accepted region {F,H,h) 


0.912 


0.667 






7r T stop without nuclear interaction 


0.492 


0.524 


MC 




or decay-m-night 










ah /ax and kinematic cuts 


0.537 


0.537 




^scat 


7r + ^ / u + ^e + cuts 


0.349 


0.349 


^ scat 




7 reconstruction and nducial cuts 


0.530 


0.492 


MC, 




7 selection cuts 


0.216 


0.177 


MC, 




Other cuts on beam and target 


0.507 


0.507 






Total acceptance 


2.99 x 10" 4 


1.10 x 10~ 4 







Table 2 

Acceptance factors for the K + — ► 7r + 77 decay in the kinematic region P > 213 
MeV/c, for c = 1.8 including unitarity corrections ("UC") and for c = 1.6 with- 
out the corrections ("w/o UC"), and the samples used to determine them. "MC" 
in the rightmost column means the sample generated by Monte Carlo simulation. 
"Kb" , "^2", "K w 2 n , and "7r sca t" mean the data samples of kaons entering the tar- 
get, K + — > decays, K n 2 decays, and scattered beam pions, respectively; these 
samples were accumulated by calibration triggers simultaneous to the collection of 
signal candidates. 
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